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Abstract: The BRST treatment of triaxial systems rotating at high spins is 
used to solve perturbatively the 7-independent Bohr collective hamiltonian. 
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O ' 1 Introduction 

G\ _ 

The standard model used to calculate nuclei at high spins is the cranking 
model . It is particularly successful in the prediction of energy levels, single- 
particle states, deformation parameters, etc. . The problem of constructing the 
normal modes associated with the cranking model has been treated by several 
authors 0. All these treatments assume the quadrupole plus monopole pairing 
interaction. A different treatment of the cranking model at high spins is 
based on the application of the BRST symmetry @|| to many-body problems 
0. Since the collective coordinates are explicitly introduced, the treatment 
requires constraints and gauge conditions. An extension of the formalism to 
include transition probabilities has recently appeared J7J. 

In the BRST treatment there are no restrictions concerning the hamiltonian[], 
other than those associated with the validity of the cranking model: small 
fluctuations in physical magnitudes relative to their expectation value (if this 
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one does not vanish). In particular, this must be true for the collective angular 
momentum components I relative to the expectation value (I x ) = I. 

None of the previous formalisms has been carried out beyond the level of con- 
structing normal modes and their corresponding transition matrix elements 
to the ground state. However (perturbative) calculations should also be per- 
formed within such basis. In the present letter we report for the first time 
on one of such calculations. It is performed by applying the formalism to a 
reasonably simple but non-trivial system, the Bohr collective hamiltonian 0. 
Our aim is: i) to verify that the problem of infrared divergencies has been 
properly solved and ii) to compare the perturbative results in the cranked 
basis with those in the laboratory system (which can be solved in this case). 
Emphasis is put in those diagrams associated with the modification of the 
rotational energies in excited states, since their systematic application may 
shed light on the still unsolved problem of identical bands. 

The Bohr collective hamiltonian is constructed out of the five components of 
a quadrupole tensor and their associated conjugate momenta p M . We prefer 
to work with the "cartesian" components q v 

«o = qp] a±i = ±\fl(q y ± iq x ) ; «±2 = \[\{<h T (1) 



In terms of such coordinates and momenta, the kinetic energy and the two 
(independent) static invariants of the problem may be written 

T=±M + P 2 1 +pI+pI+P 2 z ) (2) 

W 2 = <$ + q* + ql + q\ + q\ = [3 2 (3) 
W 3 = -2q\ + 6q p <% + 3q p (2q 2 z - q 2 x - q 2 y ) + 3V%(g* - q 2 y ) 

-QV3q x q y q z = -2(3 3 cos 3 7 (4) 



We work in a rotating frame of reference, and we formulate the problem within 
an overcomplete space consisting of the five original^ or intrinsic variables q v 
(and their conjugate momenta p v ) plus three collective coordinates, the Euler 
angles <p v (and the three components of the collective angular momentum in 
the moving frame I u , [I^Iy] = —ie^plp). The components of the intrinsic 
angular momentum are given in appendix 

Since the collective variables are raised to the level of true variables, some 
combination of the intrinsic variables belong to the spurious sector, which is 

2 These five original degrees of freedom of the Bohr hamiltonian are usually sup- 
posed to represent collective coordinates. Note that this is not the interpretation 
here, since we reserve the name collective (as opposed to intrinsic) to the Euler 
angles determining the position of the rotating frame. 
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completed through the introduction of three Lagrange multipliers Q u (and 
their conjugate momenta P u ) and six ghost (fermion) variables r\ v ,r\ v (and 
their conjugate partners 7i u ,n u ). 

The BRST hamiltonian for an ellipsoidal system reads 

1 A 

-f^BRST = Hb — £lv{Jv — If) + —G V P V — 77=7-P„ 

+VK v Tt v - rrV^G^, J„] + ie^pQp-K^TCy (5) 

where Hb is the original Bohr hamiltonian, and are gauge-fixing functions 
which we choose later. Three of the constants F u , A v disappear during the 
treatment, while the remaining three give rise to the spurious frequencies 
u) x ,u p (eqs. flr3"D and (TTS])). 

The coordinates q@, q 1 and the x-components of the intrinsic and collective 
angular momentum and of the Lagrange multipliers have the non-vanishing 
expectation values 

(qp) = A>cos7 ; (q 7 ) = /5 sin7 (6) 
(J*) = (4) = 3<n*> = 1 (7) 



The minimization of the hamiltonian -^brst together with the condition (|7|) 
implies the minimization of the routhian S((H — £IJ X )) = (O = (£l x ))> an d 
yields the value of the deformation parameters /3 (J),7 (J) and the moment 
of inertia S(J). Therefore the cranking solution is the starting point of a 
treatment of collective rotational motion which is exact in principle. 

If the potential energy surface V = V^W^W^) is independent of W3 (i.e., 
7-independent) the minimization condition yields 

7o = vr/6; 3 = A(3 2 Q ; (p x ) = 2np o (8) 

Since there is no restoring force in the 7-direction, the system adquires in- 
mediately the 7 deformation that minimizes the kinetic energyTj. 

Because W 2 = f3 2 , the value of [3 is determined from the eq. : 




3 This fact explains the difficulties in deciding empirically between the 7-unstable 
and the spheroidal rotor descriptions. 
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where Vi = d up%y \p 2 =pl- The expectation value of the BRST hamiltonian is 

<#BRST> = ^ + Wo) (10) 



In the next step we construct the normal modes. Each operator may consist 
of several terms which are labelled by the number of phonons through the 
supraindex (t) (for instance, (a) = a^). The normal modes are determined 
from the quadratic terms in the hamiltonian. Due to signature-invariance, 
these terms split into two parts which are labelled by the subindices x and _L, 
respectively. The expression for H® is 

T ax ~^^Bx ^"^ ^ xx 

= Up? 2 + p? 2 + p? 2 ) + v^T + a? 2 + <ii l)2 ) 
+\v 2 ^ q f + q ?) 2 - oj£) (ii) 

H% = -n« J?) + yG x P x - + ivr x vr x + ±- Vx f) x (12) 

H® = Hj£ + H® = tJp^Tp + §) + uj,(T% + |) 

) (13) 

The hamiltonian includes the two real vibrations with frequencies u; 7 = 2Q and 
uj/3 = V8^i + and a (supersymmetric) spurious sector arising from (12) 
and from a term (jW) 2 /^^ obtained in (|TTD . The existence of this term is a 
consequence of the fact that the first two lines commute with the operator . 
The construction of the spurious sectorQ is made as in HP,|7| . The (arbitrary) 
frequency u x should disappear from any physical result. The (linear) gauge 
G x is the RPA angle which is obtained simultaneously with the dynamical 
moment of inertia from the RPA eqs. [JT0| 



[H2,G a ] = ~jp- t [G„J«]=i (14) 




Note the presence of the quadratic term —QJW (eq. ( |A.8| )) in the hamiltonian 
with positive signature flllD , which has been overlooked in previous derivations 
of the normal modes in cranked systems. 

4 Note the relations [r ls ,rj s ] = -[r 0s ,rj s ] = [a s ,a s ] + = [b s ,b s } + = 1 for s = x,p. 



4 



The quadratic hamiltonian in the perpendicular direction is 

Hi! = W y 1)2 + p? 2 ) + vM y 1)2 + a { z 1)2 ) - n( J« - 4 2) ) (16) 

Hi! = -0,(J« - /«) + ^G S P„ - ^P, 2 - - I«) + ±G Z P Z 

y y z 

A i i 

rr(2) _ „(2) rr(2) 
-n± - #a± + H b± 

= 1^1 + uj w (TlT w + \) + a; p (rl p ri p - rj p r 0p + a p a p + b p b p ) (18) 

where p — ±. The linear and quadratic terms of the collective angular mo- 
mentum are obtained by means of the Holstein-Primakoff representation ||11|| . 
The first (constant) term in (O) plus the rotational term in eq. (ffO) yield 
the quantal form /(/ + 1). The real degree of freedom in ([18]) represents the 
wobbling motion with frequency u w = 3f2. There also appear two spurious 
sectors with (arbitrary) frequencies u p . The gauge functions and the moments 
of inertia are derived from a generalization of the RPA eqs. ([14]) |§ 

[Hau G y ] = -±- J« + il ( i- - i) G z ; [G y , J«] = i (19) 

[Hol, G z ] = ~ J« - il - ^) G y ; [G z , J«] = i (20) 

plus the condition [G y ,G 2 ] = 0. These eqs. yield G y = — Po 1( ly j = 
-/?- x g« and ^ = 3 2 = /3 2 . 

The quadratic hamiltonians ( pT3|) and (0) determine both the real and the 
spurious spectra associated with a given value of /. The corresponding states 
span a space which is factorized into real, spurious and rotational subspaces 

\np,n 1 ,n w ) I \n 0l/ ,n lv ,n au ,n bv ) I \IM) = (T^ (T\) n - (Tt) n - \r)j 

x(Tlr°" (rL)" 1 " {a v ) n ^{b v )^ \sp)j \IM) (21) 

where rip, n 7 , n w , rii u , n 0u = 0,1,2,.. ., and n au , n bv = 0, 1. The product of the 
vacua of the real and spurious excitations (|r)i and |sp)i, respectively) is also 
the yrast state for a given /. The rotational states \IM) are, in configuration 
space, the wave functions D\ n of the rigid top. The rotational excitations 
D I MK (K < I) are partly associated with the wobbling motion and partly 
with the spurious sector. Here K is the projection along the intrinsic 

We calculate the matrix elements of the quadrupole operator in the labora- 
tory frame, since they are physical operators. In order to apply the formalism 
developed in |7j it is convenient to express the intrinsic operators with respect 
to the x-axis 

«f b) = DlMDl p (-^, -ivr, OK (22) 
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If acting within states having I x m I ^> 1, the rotational matrices may be 
expressed in terms of the collective generators I v times the operators 



increasing the angular momentum of the ground state from I to I + i (12 
In such a way we can label the components of fl22D both by the number of 
phonons and by the change in /. In particular, the zero-phonon component is 

(4 lab ))(°) = -ElP /y/2 (23) 

corresponding to the transition matrix element between yrast states (which 
differ in two units of angular momentum). 

The linear transition operators can be constructed from the expressions for 
the intrinsic variables in terms of the normal modes 

(«i lab) ) (1) = ^(Tt + r 7 ) ; («f ab) ) (1) = ^ElT w 

(c4 lab) )« = -J-JLeZ (rj (up - 40) + Vp(u p + 40)) (24) 



The classification of levels in terms of the quantum numbers (rip, n w , n 7 ) allows 
for a new interpretation of the degeneracies associated with the 7-unstable 
spectrum: for each /, there is vibrational spectrum with n 1 = 0, 1, 2, . . .. After 
n 7 = 2, they become degenerate with levels with one wobbling phonon more 
than the corresponding states with I ± 1 (fig |]). 

The vanishing of the zero-phonon terms proportional to with 1 = (di- 
agonal matrix elements) and of the matrix elements increasing in one unit 
both the angular momentum and the number of wobbling bosons is due to the 



consevation of the 7-parity quantum number (see [14 



Perturbative corrections can be calculated straightfordwardly since -£^brst is 
free of zero frequency modes. The perturbative expansion yields the residual 
hamiltonian H res = -Hbrst + ^brst + ' ' "> with 

<st = «V3g« + q ^)( q f 2 + g« 2 + e )2 + 4 1)2 + ^ 1)2 ) 

+§v 3 (Vs q ? + g «) 3 - n x (4 2) - 4 2) ) - M 2) - M 2) 

-—V^r/AG^ J i 2) ] + iW°pV" ( 25 ) 

+/3 V 3 (v / 3^ ) + qUfiqW + g« 2 + + g« 2 + g« 2 ) 
+§^(v / 3g« + q^Y + 0,if + njf) (26) 



6 



The Jj) 2 ^ are given in appendix [A] and Ijp are obtained by means of the 
Holstein-Primakoff representation. 

-f^BRST an d -^brst & ve T ^ se t° the vertices shown if fig. [|. There are no vertices 
with an odd number of _L-phonons due signature-invariance, and there are no 
vertices with an odd number of ghosts due to the fact that -£^brst is an even- 
Grassman function. 

The diagrammatic corrections of o(/3 2 ) to the energies of the yrast state and 
of the one /9-phonon state are given in the second and third rows of fig. 
respectively. Although individual diagrams depend on the spurious frequen- 
cies uj x and u± this dependence cancels out from the final sums. The total 
corrections to the energies are 

AE (yrast state) = 




2 



-s^(- m + ^ + T v ^) < 27 > 

AB(one /J-pho n o„) = — + ^-V 1+ -V 2 + -V 3 + -V t j 
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The 7-independent Bohr collective hamiltonian may be analitically solved in 
the laboratory frame through the conventional method of separation of vari- 
ables ||14|| . Within such a framework (which does not require the treatment of 



a spurious sector) corrections fl27j) and ( Pq) can be reproduced analitically. 



In the case of the harmonic motion ||, V\ = constant and = 0. Conse- 
quently, the "static" moment of inertia is S = 2I/y/2V%, while the "dynam- 
ical" is $2V X — ► oo, as required to reproduce yrast energies increasing linearly 
with /. 

A similiar treatment can be applied to a 7-dependent potential, the main 
difference being that 7 Q = j (I) varies from the equilibrium value at / = 
to the value (^J) for large /. In this case the solution requires a numerical 
computation. 

We have perfomed (to our knowledge, for the first time) a perturbative cor- 
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(1,0,0) 




7-2 7-1 7 (even) 7 + 1 7 + 2 

Fig. 1. The spectrum of -^brst* ^he l eve ls are classified by (np ). Dark lines 

correspond to strong transitions between yrast states, lighter ones to transitions to 
vibrational states, and dotted lines to 7-parity inhibited transitions. 

rection to zero order results appearing through the application of the cranking 
approximation. This has been done for the Bohr collective hamiltonian with 
a 7-independent potential energy surface. Although the empirical applicabil- 
ity of such a model is being increasingly recognized, we mostly consider the 
present contribution as a non-trivial exercise for future calculations with re- 
alistic systems. In particular, the replacement of the initial and final phonons 
by a fermion line and the corresponding modification of diagrams measures 
the difference between the moments of inertia associated with the even and 
odd systems (the problem of identical bands). 



A Intrinsic angular momentum operators 

The intrinsic angular momentum operators are defined by J M = i\/l0(q, p)*. 
In cartesian components they are 

Jx = QzPy ~ %Vz + {VZqp + q^Px ~ QxiVSpp + P-y) (A.l) 

J y = q x p z - q z Px - (Vzqp - q~,)p y + q y (\^p/3 - p~,) (A.2) 
Jz = q y p x - q x p y - 2q~fPz + lq z p 1 (A.3) 
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X 




Fig. 2. The vertices and diagrammatic contributions to the energies of the yrast and 
one /3-phonon state. Here x(_L) stands for /3, 7, l x , l x (w, l p , P ) phonons and g for 
ghosts. The isolated crosses denote the expectation value of H^ ST . 



The expansion of the intrinsic coordinates and momenta gives 



(J x ) = I ; (J y ) = (J z ) = (A.4) 

4 1} = WoP? + ^(V&Z? + q?) (A.5) 

W = -Pop? - ±f ; J? = -to? + _^a> (a.6) 

^ *^xx ^ x, I * ^0 

■/S = + tfW - d^pf + P?) (A.8) 

42 = q?V? - q£¥? (A.9) 
J? = ~ q ( M ] ' (V3q? ~ tfW + - p«)(A.10) 

4 2) = 9?^ - «W " 2#W> + (A.ll) 
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